On weak vs. strong universality in the two-dimensional 
random-bond Ising ferromagnet 
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We address the issue of universality in two-dimensional disordered Ising systems, by considering 
long, finite-width strips of ferromagnetic Ising spins with randomly distributed couplings. We cal- 
culate the free energy and spin-spin correlation functions (from which averaged correlation lengths, 
£ ave , are computed) by transfer-matrix methods. An ansatz for the size-dependence of logarithmic 
corrections to £ ave is proposed. Data for both random-bond and site-diluted systems show that 
pure system behaviour (with v = 1) is recovered if these corrections are incorporated, discarding 
the weak-universality scenario. 
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It is well known that the Harris criterion (!]], for the 
relevance or irrelevance of weak disorder upon critical 
behaviour at a phase transition, is inconclusive for the 
two-dimensional Ising model where the specific heat of 
the pure system diverges logarithmically at the critical 
point. A great deal of effort has been dedicated to eluci- 
dating the properties of disordered versions of this model 
0). Both weak and strong disorder have been consid- 
ered. Early proposals implying strong deviations from 
pure system behaviour, such as a peculiar exponential 
decay of critical correlations with distance and a mag- 
netization exponent /3 = || have been ruled out by 
extensive numerical simulations 0. Though nowadays 
there seems to be general agreement that no such dras- 
tic changes are expected to arise from disorder in this 
case, two main pictures have taken hold recently, which 
seem to be mutually excludent. The first, typically rep- 
resented by the work of Heuer and of Talapov and 
Shchur P] maintains that the critical behaviour is unaf- 
fected by disorder (apart from possible logarithmic cor- 
rections, which though not explicitly considered in Refs. 
j| and H , fit in with similar overall conclusions ) ; this 
view agrees with early numerical work on magnetization 
moments || . According to the second view [Qj8| , critical 
quantities such as the zero-field susceptibility and cor- 
relation length display power-law singularities, with the 
corresponding exponents 7 and v changing continuously 
with disorder; however, this variation is such that the ra- 
tio j/v is kept constant at the pure system's value (the 
so-called weak universality scenario ||). 

The present work aims at shedding light into this 
controversy, by means of strip calculations which, since 
the work of Nightingale (l(]|ll[] connecting Finite-Size 
Scaling Jl^|l3| and Renormalization Group ideas, have 
proved to be among the most accurate techniques to ex- 
tract critical points and exponents for non-random sys- 
tems in two dimensions. Extensions of this approach to 
random systems require an appreciation of the subtleties 



involved in the corresponding averaging process pT| . [l5| ; 
early efforts in this direction [[l6) have since been ex- 
tended and put into a wider perspective Jl7],[l^,[ll| . We 
consider a two-dimensional, square-lattice, random-bond 
Ising model with a binary distribution of ferromagnetic 
interaction strengths, each occurring with equal proba- 
bility. For this specific model, the transition tempera- 
ture is exactly known from duality ||^,|2l|, so one can 
be sure that numerical errors due to imprecise knowl- 
edge of the critical point (such as may happen e.g. for 
site-diluted cases) are absent. The only sources of such 
errors will then be the finite strip width and those arising 
from the averaging process. However, the former can be 
controlled by finite-size scaling theory |ll|,0], while the 
effects of the latter are reduced by studying large enough 
samples (though this is a subtle point when correlation 
functions are concerned, as seen below). A previous strip 
calculation for this model jl9| concentrated on testing 
for random systems the well-known relation p2] ] between 
the critical exponent r\ and the correlation length on a 
strip at the critical temperature of the two-dimensional 
system, as well as on extracting the conformal anomaly 
c (proportional to the leading finite-width correction to 
the bulk free energy |^,^|). The conclusion was that, 
within error bars, r\ = 1/4 and c = 1/2 (the pure Ising 
values) for wide ranges of disorder. Since 77 = 2 — 7/f, 
those results could not be used to test any disorder de- 
pendence of 7 and v separately. Here, instead, we resort 
to numerical derivatives to obtain v. 

We have used long strips of a square lattice, of width 
4 < L < 12 sites with periodic boundary conditions. 
In order to provide samples that are sufficiently repre- 
sentative of disorder, we iterated the transfer matrix |Tl| ] 
typically along 10 7 lattice spacings, meaning much longer 
strips than those used in Ref. [jl9| . 

At each step, the respective vertical and horizontal 
bonds between first-neighbour spins i and j were drawn 
from a probability distribution 



1 



P(Jij) = ^(S(Ja - Jo) + S(Jij - rJ )), < r < 1, (1) 



which ensures POpifl that the critical temperature (3 C = 
1/ksTc of the corresponding two-dimensional system is 
given by 



sinh(2/? c J )sinh(2/? c rJ ) = 1 



(2) 



We have used three values of r in calculations: r = 0.5, 
0.25 and 0.1; the two smallest values have been chosen 
for the purpose of comparing with recent Monte-Carlo 
simulations where v and 7 are evaluated A wide 

range of disorder is thus covered. 

The procedure for evaluation of the largest Lyapunov 
exponent A" for a strip of width L and length N 3> 1 is 
well known The average free energy per site is 

then f L ave {T) = -j;A L , in units of k B T. 

From finite-size scaling, the initial susceptibility of a 
strip at the critical temperature of the corresponding in- 
finite system, xl{T c ) must vary as : 



Xl(T c ) 



d 2 f L ave (Tc) 



dh 2 



= IP> V Q(0) 



(3) 



/i=0 



where ft is a uniform external field and Q(0) is a con- 
stant [jl2| . A s f T ave (T) is expected to have a normal dis- 
tribution [ |15|j26| , so will xl- Thus the fluctuations are 
Gaussian, and relative errors must die down with sample 
size (strip length) N as 1/ \^N. The intervals (of external 
field values, in this case) used in obtaining finite differ- 
ences for the calculation of numerical derivatives must 
be strictly controlled, so as not to be an important ad- 
ditional source of errors. We have managed to minimise 
these latter effects by using Sh typically of order 10~ 4 in 
units of J when estimating f L ave (T c ; h = 0, ±8h) for the 
derivative in Eq. (||). 

A succession of estimates, (j/v) L , for the ratio 7/1/ is 
then obtained from Eq. (|^) as follows: 



In Ixl(T c )/ X l-i(Tc)} 
In [£/(£-!)] 



(4) 



Least-squares fits for plots of {^/v) L against 1/L 2 (see, 
e.g., Ref. jl9| for a discussion of suitable powers of 1/L 
for extrapolation) provide the following results: 7/1/ — 
1.748 ± 0.012, 1.749 ± 0.008, and 1.746 ± 0.013, respec- 
tively for r = 0.50, 0.25, and 0.10; the latter two esti- 
mates agree with 1.74 ± 0.03, 1.73 ± 0.05, obtained in 
Ref. || . 

The overall picture is thus consistent with 7/1/ = 7/4 
for all degrees of disorder. Taken together with the re- 
sults of Ref. jnj, and using the scaling relation 7/z/ = 
2 — 77, this confirms the view that: (1) the conformal 
invariance relation |2^] r) = L/tt^l(T c ) still holds for dis- 
ordered systems, provided that an averaged correlation 
length is used; and that (2) the appropriate correlation 



length to be used is that coming from the slope of semi- 
log plots of correlation functions against distance [ fl9| |. 

We now present results for the exponent v. The first 
difference to the free energy calculation described above 
is that the correlation functions are expected to have a 
fo<?-normal distribution Jl4||l5| rather than a normal one. 
Thus self-averaging is not present, and fluctuations for a 
given sample do not die down with increasing sample size. 
However, we have seen that overall averages (i.e. central 
estimates) from different samples do get closer to each 
other as the various samples' sizes increase. Accordingly, 
in what follows the error bars quoted arise from fluctua- 
tions among four central estimates, each obtained from a 
different impurity distribution. Similar procedures seem 
to have been followed in Monte-Carlo calculations of cor- 
relation functions in finite (L x L) systems ||. 

The direct calculation of correlation functions, (ooctr), 
follows the lines of Section 1.4 of Ref. O], with standard 
adaptations for an inhomogeneous system |]l9f . For fixed 
distances up to R = 100, and for strips with the same 
length as those used for averaging the free energy, the 
correlation functions are averaged over an ensemble of 
10 4 -10 5 different estimates to yield (croc/i)- 

The average correlation length, £ > ave , is in turn defined 

by 



(a a R ) ~ exp (-i?/£<™ e ) , 



(5) 



and is calculated from least-squares fits of straight lines 
to semi-log plots of the average correlation function as a 
function of distance, in the range 10 < R < 100. 

We can then apply the usual finite-size scaling (FSS) 
arguments |l2],[l3| to obtain estimates Vl of the exponent 
v. Assuming a simple power-law divergence - i.e., ig- 
noring, for the time being, less-divergent terms such as 
power-law or logarithmic corrections - of the correlation 
length in the form £ ~ t^, with t being some reduced 
distance to the critical point, its FSS ansatz becomes 



er = l Hz), 



(6) 



where z = tL 1 ^ and T is a scaling function. Since v does 
not appear explicitly in the expression for ££" e (T c )/L, 
one resorts to the temperature derivative of the correla- 
tion length, which can also be cast in a similar scaling 
form, 



Mi 



dt 



= L 1+ i G(z) 



with Q = dJ- ' jdz. fiL at T c (see Eq. |2|) is calculated 
numerically from values of £™ e evaluated at T c ± ST, 
with ST/T C — 10~ 3 . For systems of sizes L and L — 1, 
one obtains the estimates 



In (pl/hl-i) t= 



T=T a 



ln(L/L - 1) 



(8) 
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Note that this is slightly different from the usual 
fixed-point calculation 1 10 ll[], and is more convenient 
in the present case where the exact critical temperature 
is known. Our data for each pair of (L,L — 1) strips 
are shown in Table I, together with results of extrap- 
olations against l/L 2 for each separate sequence corre- 
sponding to different values of r. Taken at face value, 
the data show a systematic trend towards values of v 
slightly larger than the pure-system value of 1, though 
the variation is smaller than that shown in Ref. [p5| . 

Before accepting this trend as an indication of the 
weak-universality scenario, we must test for corrections 
caused by less- divergent terms as being responsible for 
the observed disorder dependence of v. We first recall 
that logarithmic corrections have already been proposed 
for the bulk correlation length in the form p] 



£~i-"[i+cki(i/t)r, 



(9) 



with v = 1 and v = 1/2, and C is a (disorder-dependent) 
constant. For the same reasons as above, estimates for v 
can only be consistently tested through the temperature 
derivative of £ : 



dt 1 



Cln(l/t)f. 



(10) 



plus less-divergent terms. 

For finite systems, logarithmic corrections are expected 
to show on scales larger than a disorder-dependent char- 
acteristic length Lq ~ exp(l/C) ||. A finite-size scaling 
ansatz for the behaviour at T c can be obtained by a suit- 
able generalisation of the standard procedures for pure 
power-law singularities (see, e.g., Ref. |^3| ). Assuming 
v = 1 and 9 = 1/2 one then has, to dominant order: 



L 2 



(l-61n£) 1/2 , 



(11) 



where b ~ 1/lnic- Figure [l] shows the results for 
(/i^/L 2 ) as a function of InL, for different values of 
r. In each case, log-corrected behaviour sets in for suit- 
ably large L, exactly in the manner predicted by theory: 
the data stabilize onto a straight line only for L > Lq, 
which decreases with increasing disorder H . 

This crossover effect is similar to that found by Wang et 
al. [|tJ in their fitting of specific heat data to the double 
logarithmic divergence predicted by theory [9. However, 
the increasing broadening of the specific heat maximum 
with disorder in finite-sized systems has been interpreted 
as evidence against double-logarithmic corrections 0,^8| . 
Here, instead, we deal with a case of single-log corrections 
to a divergence much stronger than that of the specific 
heat. It is thus easier to separate between corrections and 
the dominant power-law behaviour, as made evident by 
the consistent fits displayed in Fig. [I]. Though £,l{T) and 
the susceptibility Xl(T) were calculated through Monte 
Carlo simulations in Refs. [Mps!, no attempt seems to 



have been made to fit the corresponding data to a form 
similar to Eq. (|ll]) . It must be recalled that by examining 
the behaviour at the exact T c , only finite-size effects play 
any role in tuning the crossover; by contrast, when T c is 
not known exactly, one cannot be sure whether the (ther- 
mal) crossover towards critical behaviour has already oc- 
curred. As a final check of our data, we have also tried 
less-divergent power-law corrections, but the fittings were 
always much poorer than those assuming logarithmic cor- 
rections. In view of these facts, the estimates provided by 
Eq. (||) should then be regarded as effective exponents, 
since strong universality still holds. 

We have also applied the ideas behind Eq. ( |TT| ) to the 
two-dimensional site-diluted Ising model, using the cal- 
culational scheme proposed in Ref. (lj]. Results for the 
pure and diluted cases (for concentrations of magnetic 
sites in the range p — 0.65 - 0.95) and L = 3 — 7 are 
depicted in Fig. ||. It can be seen that the qualitative 
trend clearly changes towards a In L-dependence similar 
to that found in the random-bond case as soon as dilu- 
tion is introduced. The small curvature in the plots must 
be at least partly attributed to imprecise knowledge of 
the exact critical line (and to the approximate nature 
of that calculational scheme itself, which is asymptoti- 
cally exact only as T — > p!s[). Again, pure-system 
exponents with logarithmic corrections (as opposed to 
dilution-dependent ones |0,||) seem to describe the be- 
haviour of site-diluted Ising magnets in two dimensions. 

In conclusion, our data independently confirm that the 
conformal invariance result £ a " e = L/irrj is still valid 
for the two-dimensional random-bond Ising model, with 
r/ = 1/4 as in the pure case. The apparent dependence of 
v with disorder was found to be due to logarithmic cor- 
rections, which become more important the farther one 
moves away from the pure (i.e, r = 1) system. The weak- 
universality scenario, though quite appealing for the pos- 
sibility of demanding new underlying concepts to be ex- 
plained, docs not seem to hold in the two-dimensional 
random-bond Ising model. A similar picture most likely 
holds for the site-diluted model as well. The results pre- 
sented here do not necessarily imply, however, that con- 
formal invariance or strong universality should be valid 
for any type of disorder. In the problems treated here 
correlations can still freely propagate, unlike cases where 
e.g. frustration is allowed; we are currently investigating 
these issues for spin-glass-like systems. 
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Extrap. 
Ref. M 



1.000 ±0.030 
1.009 ±0.033 
1.012 ±0.032 
1.016 ±0.039 
1.037± 0.016 



1.052 ±0.026 

1.063 ±0.016 
1.062 ±0.017 

1.064 ±0.032 
1.083 ±0.014 

1.09 ±0.01 



1.13 ±0.06 
1.15 ±0.07 
1.15 ±0.08 

1.14 ±0.11 
1.14 ±0.06 
1.23 ±0.02 



FIG. 1. Finite-size scaling plots of logarithmic corrections 
[Eq. (111])]. Straight lines are least-squares fits of data re- 
spectively for L = 9 - 12 (r = 0.5); 7 - 12 (r = 0.25) and 
4 — 12 (r — 0.1). For data in this figure, fiL = d^/dK, with 
K = J IT. 



FIG. 2. Finite-size scaling plots of logarithmic corrections 
[Eq. ([n])] for site-diluted Ising model. Top to bottom: con- 
centration of magnetic sites = 1.0, 0.95, 0.90, 0.80, 0.70. 0.65. 
See Ref. (ll| for details on the calculation of £l. 



TABLE I. Critical exponent v from Eq. 



L r = 0.50 

5 0.928 ±0.004 

6 0.962 ±0.010 

7 0.981 ±0.020 

8 0.997 ±0.025 



0.25 0.10 

0.993 ±0.015 1.13 ±0.11 

1.029 ± 0.024 1.15 ±0.06 

1.040 ±0.036 1.14 ±0.06 

1.053 ± 0.040 1.15 ±0.06 
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